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Abstract 



After recalling the precise existence conditions of the zeta function of a pseudodifferential 
operator, and the concept of reflection formula, an exponentially convergent expression for 
the analytic continuation of a multidimensional inhomogeneous Epstein-type zeta function 
of the general form 

HeZ p 

with A the p x p matrix of a quadratic form, bap vector and q a constant, is obtained. It 
is valid on the whole complex s-plane, is exponentially convergent and provides the residua 
at the poles explicitly. It reduces to the famous formula of Chowla and Selberg in the 
particular case p = 2, b = 0, q = 0. Some variations of the formula and physical applications 
are considered. 
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1. Introduction: Existence of the zeta function of a pseudodiffer- 
ential operator ($DO) 

A pseudodifferential operator A of order m on a manifold M n is defined through its symbol 
a(x,£), which is a function belonging to the space S m (R, n x R n ) of C°° functions such that 
for any pair of multi-indexs a,/3 there exists a constant C a ^ so that 

^a(x,0|<C^(l + |^|r-H. (1) 

The definition of A is given, in the distribution sense, by 

Af(x) = (2vr)- J e^>a(x, 0/(0 ^, (2) 

where / is a smooth function, / e S [remember that S = |/ G C 0O (R n ); sup :r |x /3 (9 a /(x)| < oo, 
Va, /3 G R"}], S' being the space of tempered distributions and / the Fourier transform of 
/. When a(x, is a polynomial in £ one gets a differential operator. In general, the order 
m can be complex. The symbol of a \I/DO has the form 

a(x,0 = a m (x,0 + a m „i(x,0 H h a^^x,^) H , (3) 

being a k (x,£) = b k (x) £ k . 

Pseudodifferential operators are useful tools, both in mathematics and in physics. They 
were crucial for the proof of the uniqueness of the Cauchy problem U and also for the proof 
of the Atiyah-Singer index formula 0. In quantum field theory they appear in any analytical 
continuation process (as complex powers of differential operators, like the Laplacian) . And 
they constitute nowadays the basic starting point of any rigorous formulation of quantum 
field theory through microlocalization, a concept that is considered to be the most important 
step towards the understanding of linear partial differential equations since the invention of 
distributions ||. 

For A a positive-definite elliptic \I/DO of positive order m e R, acting on the space of 
smooth sections of an n-dimensional vector bundle E over a closed, n-dimensional manifold 

M, the zeta function is defined as 

rt 

( A (s) = tr A- s = A7 S , Re s > - = s . (4) 

The quotient s$ = dimM/ordA is called the abscissa of convergence of (a{s), which is 
proven to have a meromorphic continuation to the whole complex plane C (regular at Sq), 
provided that the principal symbol of A (that 0) admits a spectral cut: L 9 = 



{A G C; ArgA = 9, 6>i < 6 < # 2 }, SpecAnLg = (the Agmon-Nierenberg condition). Strictly 
speaking, the definition of (a(s) depends on the position of the cut Lg, not so that of the 
determinant B det^ A = exp[— <^(0)], which only depends on the homotopy class of the cut. 
The precise structure of the analytical continuation of the zeta function is known in general 
II . The only singularities it can have are simple poles at 

Sk — (n — k)/m, k — 0, 1, 2, . . . , n — 1, n + 1, (5) 

The applications of the zeta-function definition of a determinant in physics are very impor- 
tant [171 M. A zeta function with the same meromorphic structure in the complex s-plane 
and extending the ordinary definition to operators of complex order m 6 C\Z (it is clear 
that operators of complex order do not admit spectral cuts), has been obtained in Ref. ||. 
The construction starts there from the definition of a trace, obtained as the integral over 
the manifold of the trace density of the difference between the Schwartz kernel of A and the 
Fourier transformed of a number of first homogeneous terms (in £) of the usual decomposition 
of the symbol (^j of A. 

2. Exponentially convergent analytic continuation (with explicit 
poles and residua) of our zeta function 

A fundamental property shared by zeta functions of any nature is the existence of a 
reflection formula. For a generic zeta function, Z(s), it has the form: 

Z(w-s) = F(w,s)Z(s), (6) 

and allows for its analytical continuation in a very easy way — what is the whole story of the 
zeta function regularization procedure (or at least the main part of it). But the analytically 
continued expression thus obtained is just another series, which has again a slow convergence 
behavior, of power series type (actually the same that the original series had, in its own 
domain of validity). Some years ago, the notorious mathematicians S. Chowla and A. Selberg 
found a formula, for the case p = 2, b = and q = above flT0|, that yields exponentially 



quick convergence, and not only in the reflected domain. They were extremely proud of that 
formula — as one can appreciate by reading the original paper, where actually no hint about 



its derivation was given. In Ref. |L1] we generalized this expression to inhomogeneous zeta 
functions, but staying always in two dimensions (p = 2), for this was commonly believed to 
be an unsurmountable restriction of the original formula (see, for instance, Ref. |[L2|| ). The 



generalization we carried out in [11] was already non trivial and a very detailed account of 
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this step is given in |13| (where a missprint of the original derivation is corrected). Finally, 
we have realized that an extension to an arbitrary number of dimensions is actually possible. 
In fact, here we shall obtain formulas for arbitrary p and arbitrary values of b and q. 

The starting point will be Poisson's resummation formula in p dimensions, which arises 
from the distribution identity 

£ 5(x-n)= £ e i2M . (7) 

(We shall indistinctly write m • x = m T x in what follows.) Applying this identity to the 
function 

f(x) = exp (^-hc T Ax + iFxj , (8) 
with A an invertible p x p matrix, and integrating then over x G R p , one obtains 



> exp — n An + b n] = , > exp - o + 27rim A [b + 27rim 

^- V 2 y v/delA^p ; V 



neZ p 

We are going to consider the following zeta function 

1 _ -k T 



(9) 



CW*)=£ 2^ + ^ A{n + c) + q = J2 [Q{n + c) + qr\ Re (s) > £. (10) 
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The aim is to obtain a formula that gives (the analytical continuation of) this multidimen- 
sional zeta function in terms of an exponentially convergent multiseries and which is valid 
in the whole complex plane, exhibiting the singularities (simple poles) of the meromorphic 
continuation — with the corresponding residua — explicitly. The prime in the summatories 



of Eq. (1C) means that the point n = is to be excluded from the sum. Such restriction is 
irrelevant as long as q 7^ (the contribution of this single point being immediately obtain- 
able), but to define the zeta function in this way is essential in order to be able to reach the 
limit q — > (as we shall do later). The only condition on the matrix A is that it corresponds 
to a (non negative) quadratic form, that we call Q. The vector c is arbitrary, while q will 
(for the moment) be a positive constant. 

2.1. The main expression (inhomogeneous case). Use of the Poisson resummation formula (|9]) 
yields, after some work, the following expression: 

(27r)P/y/ 2 -* T(s-p/2) 2 g / 2 +P/ 4+ Vg- g / 2+ P/ 4 



cos(27rm • c) (n^A^mj K p/2 - s ( 2tt ^/2qm T A~ l m ) , (11) 

/2 



meZ p 
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where K v is the modified Bessel function of the second kind and the subindex 1/2 in Z^ 2 
means that only half of the vectors m G Z p intervene in the sum. That is, if we take 
an m G Z p we must then exclude — rh (as simple criterion one can, for instance, select 
those vectors in Z p \{0} whose first non-zero component is positive). Eq. ([11]) fulfills all 
the requirements demanded before. It is notorious to observe how the only pole of this 
inhomogeneous Epstein zeta function appears explicitly at s — p/2, where it belongs. Its 
residue is given by the formula: 

(2tt) p/2 

With a bit of care, it is relatively simple to obtain the limit of expression (|TT|) as q — > 0. 

However, instead of proceeding in this way (what we shall do later), it is adviceable to 
construct first a direct recurrent formula for the case q = 0. This is certainly the natural 
option in such case, where no cut-off q exists to safeguard the t-integration (there is no way 
to use the Poisson formula on all p indices of n at once). However, we can still deal with this 
case by using the Poisson resummation formula on some of the p indices n only, say on just 
one of them, n\. Poisson's formula on one index reduces to the celebrated Jacobi identity 
for the #3 function 



6 3 (z\t) = 1 + 2]T<f 2 cos(2nz), q = e wiT , \q\ < 1, r G C, (13) 

n=l 

the identity being: 

yj—ir 



-1 

T 



(14) 



OO 1 oo 

V"^ ^n 2 iriT+2niz ^(z—nir) 2 /(TriT) (15) 

i=— oo V n=— oo 



or, in other words 



(for a classical reference see, e.g. Ref. |1J])- Here z and r are arbitrary complex, z,r G C, 
with the only restriction that Im r > (in order that \q\ < 1). For the applications, it turns 
out to be better to recast the Jacobi identity as follows (with nir — > —t and z — ► nz): 

-n 2 t+2mnz _ \ " -ir 2 (n-z) 2 /t 



n=— 00 



equivalent ly 



E 



e 



-(n+z) 2 t 



t 

with z,t G C, Re t > 0. 



1 + >T e"" 2 " 2/t cos(2vr^) 



n=l 



(17) 
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2.2. Recurrent expression (and the homogeneous case). Using this last formula on the first 
component, n\, of the summation vector n, we obtain (for the sake of simplicity of the final 
expressions, we shall just consider now the case c = 0, but the result can be generalized to 
quite easily) : 



Ca,( s ) = 2 E H + q ) ' + th ^ v r(s " 1/2) ( H 2\-A + q 



n - ! ^\ S J _ y 



1/2- s 



2^ cos I b n 2 \n l [n 2 \-in 2 + q) K s _ 1/2 y/r4 A p _in 2 + g 



aS /2+l/4 ^ y 
ni=l 

which can be written as 

fWT(s- 1/2) 47r " NT 

Ca,0,^ S ) ~ C,0,q( S ) + y - CA p _ x ,Q,g( S - 1 / 2 ) + aS /2+l/4 r ( s ) 2-r ^ 

2^ cos b n 2 j n x \n 2 i\ p -in 2 + qj A s _i/ 2 ^— J^y n 2 A p -in 2 + qj ■ (19) 

This is clearly a recurrent formula in p, the number of dimensions, the first term of the 
recurrence being 

+^-a- 1 / 4 -/2 g i/4-/ 2 g n -V^ s _ (2™*/^) . (20) 
r(s) V V a7 

To take in these expressions the limit g — > is immediate. One obtains: 



Z P-l 



cos \^—b n 2 J n x (n 2 A p _ 1 n 2 J A s _i/ 2 ^— ^=-Y n 2 A p-i n 2j ■ (21) 

In the above formulas, A is a p x p symmetric matrix A = (tty)ij=i2 » = the 
(p — 1) X (p — 1) reduced matrix A p _i = {a i j) i . =2 , a the component a = an, 6 the p — 1 
vector b = (a 2 \, . . . , a p \) T = (a± 2 , . . . , ai p ) T , and finally, A p _i is the following (p — 1) x (p— 1) 
matrix A p _i = A v _\ — -^b®b- 

Let us now get back to the case q = starting from the beginning (e.g., the zeta function 
given by Eq. (|10|), with q — 0). From that expression, the recurrence ( ^IP can be obtained 
directly in the same way — this is rather obvious, since q plays no role in the derivation. 
What is not so obvious is to realize that the limit as q — > of Eq. ( ITT] ) is exactly the recurrent 
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formula (pip . More precisely, what is obtained in the limit is the reflected formula which 
one gets after using the well known Epstein zeta function reflection 



71-2S-P/2 



T(s)Z(s; A) = -j==T(p/2 - s)Z(p/2 - s; A' 1 ), (22) 



being Z(s; A) the Epstein zeta function |]15 |. After some thinking, such result is easy to 



understand. Summing up, we have thus checked that our formula fllTD is valid for any q > 0, 
since it contains in a hidden way, for q = 0, the recurrent expression ([21]). 

As announced at the beginning, the formulas derived here can be considered as gen- 
eralizations (in more than one sense) of the Chowla-Selberg (CS) formula. All share the 
same properties that are so much appreciated by number-theoretists as pertaining to the 
CS formula. In a way, these expressions can be viewed as improved reflection formulas for 
zeta functions; they are in fact much better than those in several aspects. Namely, while a 
reflection formula connects one region of the complex plane with a complementary region 
(with some intersection) by analytical continuation, the CS formula and our formulas are 
valid on the whole complex plane, exhibiting the poles of the zeta function and the corre- 
sponding residua explicitly. Even more important, while a reflection formula is intended to 
replace the initial expression of the zeta function — a power series whose convergence can be 
extremely slow — by another power series with the same type of convergence, it turns out 
that the expressions here considered give the meromorphic extension of the zeta function, 
on the whole complex s-plane, in terms of an exponentially decreasing power series (as was 
the case with the CS formula, that one being its most precious property). 

Actually, exponential convergence strictly holds under the condition that q > 0. However, 
the formulas themselves are valid for q < or even complex. What is not guaranteed for 
general q GC is the exponential convergence of the series, nor its power- like convergence, 
for that matter. Those analytical continuations in q must be dealt with specifically, case by 
case. The physical example of a field theory with a chemical potential falls clearly into this 
class. 

2.3. Particular case p = 2. The above statements apply both for the cases q > and q = 0. 
The last situation is more involved, however. One is led to employ the recurrence relation 
( |2T|) several times (p — 1, in general), what gives rise each time to an additional series of 
Bessel functions K v (exponential convergence). For completeness, let us write down the 



corresponding series when p = 2 explicitly. They are, with q > | T3[ 

2nq 1 ~ s 4 
Ce(s; a, b, c; q) = -q~ s + ~ _ + 



1)VA T(s) 



2^n I K s _ 1/2 2vmW- 



qa I — , V v a 



(7 



+ ^/-(2-rE^ 1/2 cos(™Va)E^ 2s (A + ^) 1/4 " S/ V s _ 1/2 ^Va + ^ 



fl n=l d|n 



where A = 4ac — 6 2 > 0, and, with q = 0, the CS formula flCfj 



s-1 



a, 6, c; 0) = 2((2s) a~ s + T ^ )As . 1/2 r(s - l/2)C(2s - 1) 

os+5/2 s oo . s 

+ r(g)A ,/2- 1 /4 v ^ E «- 1/a ^(«) cos(7rn6/a) K s _ 1/2 (_Va) . (24) 
where cr s (n) = X^|n<^ S ; sum over t ne s-powers of the divisors of n. (There is a missprint in 



the transcription of formula (|24|) in Ref. [fL2j|). We observe that the rhs's of (^) and ( |24|) 
exhibit a simple pole at s = 1, with common residue: 

Res s= iCi?(s; a, b, c; q) = —= = Res s= iCi?(s; a, b, c; 0). (25) 

V A 



3. The case of a truncated range 

The most involved case in the family of Epstein-like zeta functions corresponds to having 
to deal with a truncated range. This comes about when one imposes boundary conditions 



of the usual Dirichlet or Neumann type fL3[ . Jacobi's theta function identity and Poisson's 
summation formula are then useless and no expression in terms of a convergent series for the 
analytical continuation to values of Re s below the abscissa of convergence can be obtained. 
The best one gets is an asymptotic series expression. However, the issue of extending the 
CS formula or, better still, the most general expression we have obtained before for inhomo- 
geneous Epstein zeta functions, is not an easy one. This problem has seldom (if ever) been 
properly addressed in the literature. 

3.1. Example 1. To illustrate the issue, let us consider the following simple example in one 
dimension: 

oo _ 

( G (s;a,c;q)= £ [a(n + cf + q] \ Res > 1/2. (26) 

n=— oo 

Associated with this zeta functions, but considerably more difficult to treat, is the truncated 
series, with indices running from to oo 

oo _ 

t Gt (s-a,c-q) = £ [a(n + c) 2 + q]~ S , Res > 1/2. (27) 

n=0 

In this case the Jacobi identity is of no use. How to proceed then? The only way is to 
employ specific techniques of analytic continuation of zeta functions The usual method 



S 



involves three steps [T1J. The first step is elementary: to write the initial series as a Mellin 
transformed one 

OO _ 1 OO .QQ 

J2[a(n + c) 2 + q] dtt^expl-lain + cf + qjt}. (28) 

n=0 1 \ S ) n=0 J ° 

The second is to expand in power series part of the exponential, while leaving always a 
converging exponential factor, 

oo i oo f(XJ oo l_ n \m 

The third and most difficult step is to interchange the order of the two summations — with 
the aim to obtain a series of zeta functions — what means transforming the second series 
into an integral along a path on the complex plane, that has to be closed into a circuit 
(the sum over poles inside reproduces the original series), with a part of it being sent to 
infinity. Usually, after interchanging the first series and the integral, there is a contribution 
of this part of the circuit at infinity, what provides in the end an additional contribution to 
the trivial commutation. More important, what one obtains in general through this process 



is not a convergent series of zeta functions, but an asymptotic series [fL3[]. That is, in our 
example, 

(-a)"T(m + s) 

\ayu^c) -t-t/j ~ ^ 

n=0 m=0 

Being more precise, as outcome of the whole process we obtain the following result for the 
analytic continuation of the zeta function |L7| : 



£ \a(n + c) 2 + q] fw m+s Cg(-2m, c) + additional terms. (30) 

^ i J ^ m\T(s)q m+s 



r( 9rrY 2) g 1/2 "' + ^- 1/A ~ S/ V /A - S/2 £ cos(2 7 rnc)^_ 1/2 (2 7 ™ v ^). 
V a 2r(s) T(s) ^ v 

(Note that this expression reduces to Eq. (|20|) in the limit c — >• 0.) The first series on the 
rhs is asymptotic J16, Observe, on the other hand, the singularity structure of this zeta 
function. Apart from the pole at s — 1/2, there is a whole sequence of poles at the negative 
real axis, for s = —1/2, —3/2, —5/2, . . ., with residua: 

(2j — 1)!! q j 

Res s=1 /2-j(G t (s; a, c; q) = — — " , j = 0,1,2,... (32) 
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3.2. Example 2. As a second example, in order to obtain the analytic continuation to Re 
s < 1 of the truncated inhomogeneous Epstein zeta function in two dimensions, 



Ce t (s; a, b, c; q) = V] (am 2 + bmn + cn 2 + q)~ 

m,n=0 



(33) 



we can proceed in two ways: either by direct calculation following the three steps as explained 
above or else by using the final formula for the Epstein zeta function in one dimension 
(example 1) recurrently. In both cases the end result is the same: 



Cs t (s; a, b, c; q) = (am 2 + bmn + cn 2 + q) 

m,n=0 

(Aa) s ^ (-l)"T(m + s 



T(s) 



E - 

m,n=l 



bq 



(s-i)Ar(s- i) 



mi 

oo 



(2a)' m (A^ + 4ag)- m " s CH -2m: 



on" 
2a 



;-l) n r(ra + s - l)B n fAaq\~ n q~ s 



n=0 



nl 



V A J 



+ 



nq 



7T 

a 

1/4 



r( a - 1/2) 1/a _. 



7T 



ga 



4 2(s-l)VA 

OO y 

E ^ S " 1/2 ^-i/2 (2 



7T, 



gA 



n=l 



A 



n=l 

a 



lTm\ — 



oo 



+ W±(2tt) s ^ n s - 1/2 cos(7m6/a) £ d 1 ^ 5 A + 



4ag\ 



l/4-s/2 



?A / 



n=l 
7T7i 

a \ 



/ag 



A + 



4ag 



A 

(34) 



* " n=l d| n 

The first series on the rhs is in general asymptotic, although it converges for a wide range of 
values of the parameters. The second series is always asymptotic and its first term contributes 
to the pole at s = 1. As in the case of Eq. fl23|), the pole structure is here explicit, although 
much more elaborate. Apart from the pole at s = 1, whose residue is 

ii b 



Res s= i( Et (s;a,b } c; q) 



A' 

there is here also a sequence of poles at s = ±1/2, —3/2, —5/2, . . ., with residua: 



Res s=1 / 2 -j(Et(s; a, 6, c; q) 



(2j-l)\\qi 



1 1 

+ 



j = 0,l,2,... 



(35) 



(36) 



j!2i+ 2 Vv 7 ^ V~c) 
The formula above, Eq. (|34"D , is really imposing and hints already towards the conclusion 
that the derivation of a general expression in p dimensions for the zeta function considered 
in Sect. 2 but with a truncated range is not an easy task. 

4. Some uses of the formulas 

These formulas are very powerful expressions in order to determine the analytic structure 
of generalized inhomogeneous Epstein type zeta functions, to obtain specific values of these 
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zeta functions at different points, and from there, in particular, the Casimir effect and heat 
kernel coefficients, and also in order to calculate derivatives of the zeta function, and from 
them, in particular, the associated determinant. Notice that obtaining derivatives of the 
formulas in Sect. 2 presents no problem. Only for truncated zeta functions (Sect. 3) the 
usual care must be taken when dealing with asymptotical expansions. We shall illustrate 
these uses with three specific applications. 



4.1. Application 1. In a recent paper by R. Bousso and S. Hawking [[19], where the trace 
anomaly of a dilaton coupled scalar in two dimensions is calculated, the zeta function method 
is employed for obtaining the one-loop effective action, W, which is given by the well known 
expression 



Ci(0)lV + CA(0)J, (37) 
with Ca( s ) = t r Ar s . In conformal field theory and in a Euclidean background manifold of 



toroidal topology, the eigenvalues of A are found perturbatively (see [0), what leads one 
to consider the following zeta function: 

oo 

Ca(s)= £ (Ah)" s , (38) 

k,l=— oo 

with the eigenvalues A k i being given by 

2 2 

A w = k 2 + I 2 + — + —4 -, (39) 

2 2(4/ 2 -l)' v ; 

where e is a perturbation parameter. It can be shown that the integral of the trace anomaly 
is given by the value of the zeta function at s = 0. One barely needs to follow the several 
pages long discussion in [|T0f| , leading to the calculation of this value, in order to appreciate 



the power of the formulae of the preceding section. In fact, to begin with, no mass term 
needs to be introduced to arrive at the result and no limit mass — > needs to be taken later. 
Using binomial expansion (the same as in Ref. |19|), one gets 

oo / 2\ ~ s 2 oo / 2\ ~ 1~ 6 

C{s)= £ r + /2 + Y ~T E ( k2 + l2 + ^) (4* 2 -!)- 1 - (40) 

k,l=—oo \ / k,l=—oo \ / 

From Eq. (^) above, the first zeta function gives, at s = 0, exactly: —7te 2 /2. And this is 
the whole result (which does coincide with the one obtained in [0), since the second term 
has no pole at s = and provides no contribution. 

4.2. Application 2. Another direct application is the calculation of the Casimir energy density 
corresponding to a massive scalar field on a general, p dimensional toroidal manifold (see 

11 



||20|| ). In the spacetime M = R xE, with £ = [0, l] p /~, which is topologically equivalent 
to the p torus, the Casimir energy density for a massive scalar field is given directly by Eq. 



(0) a t 



-1/2, with q 



m 



(mass of the field), 6 = 0, and A being the matrix of the 



metric g on S, the general p-torus: 



E 



c 



•1/2). 



(41) 



The components of g are, in fact, the coefficients of the different terms of the Laplacian, 
which is the relevant operator in the Klein-Gordon field equation. The massless case is also 
obtained, with the same specifications, from the corresponding formula Eq. (j2i~P . In both 
cases no extra calculation needs to be done, and the physical results follows from a mere iden- 
tification of the components of the matrix A with those of the metric tensor of the manifold 
in question |p0[] . Very much related with this application but more involved and ambitious 
is the calculation of vacuum energy densities corresponding to spherical configurations and 
the bag model (see [21], |22], and the many references therein). 

4.3. Application 3. A third application consists in calculating the determinant of a differential 
operator, say the Laplacian on a general p-dimensional torus. A very important problem 
related with this issue is that of the associated anomaly (called the multiplicative or non- 
commutative anomaly) |[23|| . To this end the derivative of the zeta function at s = has to 
be obtained. From Eq. (|TT[) , we get 



4(2g)f/ 4 



E 



cos(27rm • c) 



VdetA ^ v (m T A- l m) 



27rJ2qrh T A~ 1 m 



(2ny/ 2 T(-p/2)qP/ 2 



+ 



VdetA 
-l) k (2Tr) k q k 



p odd, 



[ty(k + 1) + 7 — Inq] , p = 2k even, 



(42) 



fc! VdetA 

and, from here, det A = exp —(^(0). For p = 2, we have explicitly: 



det A(a, 6, c; q) 



3 27r(g-ln q)/V~A 



-27Ta/ q/a 



exp I -4 ~ 



71=1 



n 



+ cos(irnb/a) ^^d exp ( — 

d\n 



irn 
a 



Aaq 



(43) 



In the homogeneous case (CS formula) we obtain for the determinant: 



det A (a, 6, c) = - exp 
a 



-4C'(0) 



7T 



6a 



n=l 



n 



(44) 
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or, in terms of the Teichmiiller coefficients, T\ and T2, of the metric tensor (for the metric, 
A, corresponding to the general torus in two dimensions): 



detA(n,r 2 ) = ^^exp 



(45) 



Needless to mention, all the good properties of the expression for the zeta function are just 
transferred to the associated determinants, which are thus given, on its turn, in terms of 
very quickly convergent series. 
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